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Abstract 

We study parametric level statistics of the discretized excitation spectra 
inside a moving vortex core in layered superconductors with impurities. 
The universal conductivity is evaluated numerically for the various values 
of rescaled vortex velocities k from the clean case to the dirty limit case. The 
random matrix theoretical prediction is verified numerically in the large k 
regime. On the contrary in the low velocity regime, we observe Gxx oc k?/^ 
which is consistent with the theoretical result for the super-clean case, where 
the energy dissipation is due to the Landau-Zener transition which takes place 
at the points called "avoided crossing" . 
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I. INTRODUCTION 



The discrete energy spectra inside the vortex core in two-dimensional superconductors 
play an important role for transport properties in the mixed state especially in low tem- 
peratures. The energy dissipation due to the vortex motion results from the transition 
between these levels. The I-V characteristic in superconductors has been studied by many 
authors. In the superclean case ujqt 3> 1, where ujq is the level separation and is the 
inverse scattering time, the dissipative conductivity has been obtained quasiclassically as 

— eng/ {Buqt). fl]] For the 2D superconductor in the superclean limit where there is only 
one impurity inside the vortex core, the energy dissipation due to the Landau-Zener transi- 
tion between levels has been evaluated. The energy transfer takes place when the impurity 
passes through the "dangerous region" , where many Landau-Zener transitions contributing 
to the dissipation. There is also a perturbational approach for the moving impurity which 
is valid in the high velocity limit v ^ vf{^/ Ep)"^ where v is the vortex velocity and A is 
the mean level spacing. The transition between time-dependent levels Ei{X), where X is 
the external perturbation (the position of impurities) are treated in ref. in the framework 
of random matrix theory, and almost the same result for axx as the quasi-classical one is 
obtained. 

The parametric level statistics for the energy spectra under various external perturbations 
provide a new perspective of transport properties in disordered quantum systems such as 
metallic grains or quantum dots. The Thouless formula which relates the conductance to an 
equilibrium quantity such as level curvatures gave the first clue to investigations of universal 
behaviors of level statistics. P,^] The random matrix theory (RMT) and its classification 
of distributions of level spacing (Wigner-Dyson statistics), which plays a role of a finger- 
print also for quantum chaotic systems, supplies a powerful tool for the study of disordered 
systems in the framework of level dynamics. []^,|^ 

There have been both theoretical and numerical studies within the framework of the RMT 
for the level dynamics of the disordered systems such as an ring subject to the Aharonov- 
Bohm flux and non-interacting electrons in background potentials. @-12| For the quantum 
chaotic systems, there are experimental studies of the microwave billiards where the Gaussian 
distribution of level velocities has been observed. [|T3| The motivation of the current study is 
to investigate the universal behavior of the level velocities for low energy excitation spectra 
inside a moving vortex core dragged through a layered superconductor over a wide region 
from moderately clean to dirty phase. In the presence of disorder, the Hamiltonian for the 
scattering quasi-particles between the excitation levels inside the vortex core becomes the Lie 
algebra Sp{2N), and it belongs to the class C of the classification of the Atland-Zirnbauer. 
The level statistics over the crossover region from the moderately clean to the dirty 



phase has been investigated in the framework of the RMT |]T5|, and also verified numerically 



for various impurity concentration inside the core. For the energy levels far from the Fermi 
energy Ei ^ Ep in the dirty limit case, the level statistics have been shown to be identical 
with that of the Gaussian unitary ensemble (GUE). 

In section II, we investigate the universal distribution of the level velocities for the quasi- 
particle spectra within the moving vortex core for different regimes which are characterized 
by the parameter uot (or the number of impurities within the core Ni), with constant vor- 
tex velocity v. We observe that in the case of UoT -C 1 the Gaussian distribution of level 
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velocities is observed. In the clean case, non-Gaussian distribution with an enhanced peak 
at V = which is the characteristic of the systems with localized wave-functions |12| is 
obtained. 

In the diffusive time-evolution of energy spectrum within the core, the vortex velocity v 
plays an essential role. It should be large enough that pinning is negligible, since v is taken 
as a constant. However it is assumed to be moderately small f <^ f i^- in order to satisfy the 
adiabaticity of transitions between levels. The parameter vk is determined by the elastic 
mean free path and it manifests the mesoscopic nature of the level dynamics. 

In section III, the energy dissipation for the low velocity regime is investigated and 



nonlinear I-V characteristics are studied. We follow the same procedure as in ref. |T7 
in which the energy dissipation has been evaluated through iterative calculations of time- 
dependent spectra from the Schrodinger equation for a perturbed Hamiltonian. Our result 
coincides with that obtained for the superclean case in ref. in the small v regime. In the 
dirty phase {uiqT -C 1), our numerical result is identical to that of the random matrix theory 
for the Gaussian unitary ensemble with moderately large v. 

II. UNIVERSALITY IN LEVEL VELOCITY DISTRIBUTION 

The excitation spectrum inside the vortex core is obtained from the Bogoliubov-de 
Gennes equation, 



(2.1) 



where 

^ _(v^l2m-Ep^V{r) A(r) \ 

\ A*(r) -v^l2m^Ep-V[r))- ^ ' 

We consider the case for the s-wave gap parameter 

A(r) = A(r)e*^, (2.3) 

where (/) is the polar angle. With the absence of impurity {V{r) = 0), the excitation energy 
spectrum 

E^^ = -{n-^)huJo, n = 0,±l,±2,--- (2.4) 
is obtained. We set ^ = 1 hereafter and 

/ rfre-2^("), (2.5) 

JO 



K{r) = I dr'A{r')/vF. (2.6) 



r 







If we take A(r) = A (= Aqo) then the level spacing becomes cuq ~ /\ /Bp- The n-th 
eigenfunction is 
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■'n 



which we call Caroli-de Gennes-Matricon (CDM) basis where Jn{x) denotes a Bessel function 
and c is a normalization constant. 

When the impurity scattering is present, the wave function becomes a linear combination 
of the CDM basis as 

V'n = X! fmn<f>m (2.8) 
m 

where n-th column of the matrix / is given by the first order perturbation as 

ii-E + J2^')fn = 0, (2.9) 

i 

^nm ^nmE^j (2.10) 

Km = l".e-2^('^>)+*'^»('-")[J„(fcpr,)J„(A;pr,) - Jn-i{kFr^)Jm-i{kFr^)]/\F^, (2.11) 

with the assumption that impurity potential V{r) to be a sum of short-range scattering sites 
ri = {ri, (pi) 

V{r) = Y.V,6{r-ri). (2.12) 

i 

The summation for A* in eq. (|2.9| ) is taken for all impurity sites and ^ is the coherence length. 
When the vortex is dragged within the sample due to the Lorenz force, the Hamiltonian 

n{x{t)) = i + J2^' (2.13) 

i 

is perturbed by changing realizations of the impurity positions within the core. The per- 
turbation parameter X{t) is the vortex coordinate and the energy levels go up and down 
at different times (fig.l). Then the Hamiltonian becomes time dependent with the vortex 
velocity v fixed as a constant. 

The universal Gaussian distribution of the level velocities was studied for disordered 



metallic rings subject to the Aharonov-Bohm flux. [|TT| The level dynamics become universal 
dependent only on the Wigner-Dyson ensemble after a rescaling of the energy E and the 
controlling parameter X as 

six) = E{X)/A, (2.14) 
X = ,JC{0)X, (2.15) 

where A is the mean level spacing and C(0) is the normalized dispersion of level velocities. 
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The statistical average < • ■ ■ > is taken over samples and X, i.e. over time steps as we fix 
the vortex velocity v. It was shown that the equilibrium quantity C(0) is directly related to 
the transport (dissipative) conductivity given by Kubo formula. Thus C(0) is also called 
the universal conductance. 

The density of states of quasi-particles inside the vortex core in dirty or moderately clean 
phase has been investigated theoretically. |l^j2^|T^ It was also shown from the explicit 



diagonalization of the Hamiltonian matrix that in the dirty limit the Wigner-Dyson statistics 
of level distributions become identical to that of the Gaussian unitary ensemble when the 
levels are far from the Fermi energy. |T^| Following the same procedure we evaluate the 



distribution of the rescaled level velocities P{ei) where ii is written as 

dE, 1 



dt 



AJCiO) 



(2.17) 



for various numbers of impurities Ni within the vortex core. 
The inverse scattering time is given by 

r-^ = 2niV-^m. (2.18) 

rii = Ni/S,'^ is the concentration of the impurities within the core of radius S,-, where ^ is the 
coherence length. In the numerical calculations we set = 1 and ^kp = 50. Thus from the 
relation C, ~ kp/mA we have 

iuoTy' = Ni. (2.19) 

The size of the Hamiltonian matrix 7i is taken as 100 x 100. The vortex motion is restricted 
in y-direction and periodicity in this direction for the configuration of moving impurities is 
assumed. The discretized shift Ay of an impurity at discretized time step At is Ay = ^vAt. 
The rescaled level velocity ii for 100 time steps are calculated and the P{ei) for levels between 
EiQ to E40 are averaged, where we denote Eq for the Fermi Level [E^ = Ep). The sample 
averaging has been taken over 320 (480 for Ni = 500 case) different initial configurations of 
the impurities. 

In fig. 2 we show averaged distribution of the level velocities P{£) for various numbers 
of impurities within the core. We set v = 0.002 and At = 1. The inset shows the rescaled 
variance C(0) for each energy level. The condition for the dirty phase 1 ^ {Xp/O'^NNi <^ 
where N is the half of the size of the Hamiltonian is well satisfied in the case Ni = 500 
{uoT = 0.002). In this dirty phase the Hamiltonian can be described as the random matrix 
and the Gaussian distribution of the level velocities are observed. 

In the clean case Ni = 3 {ujqt ^ 0.33), the non-Gaussian behavior is found. We obtain a 
sharp peak at £ = which is similar to the result of the microwave experiment for the case 
when the local level dynamics are concerned. W^ This non-Gaussian behavior suggests 



that the wave functions for the scattering quasiparticles in this regime is still localized. |12 



III. DISSIPATIVE PROPERTY OF THE WAVE FUNCTION OF 
QUASI-PARTICLES INSIDE A VORTEX CORE 

We evaluate the evolution of the wave function \ipit)) = {\ipiit)) , \ilj2(t)) , ■ ■ ■ , \'^2Nit))} 
from the time dependent Schrodinger equation 
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ij^\m) = n\m), (3.1) 

where 

H = i + J2^' (3.2) 

i 

or equivalently consider time evolution of the coefficient matrix / which depends on time. 

= nf{t) (3.3) 

The initial state /(O) is assumed to be the i-th eigenvector of Ti. which corresponds to the 
i-th energy level Ei. 

/(0) = {0,...,/„...,0} (3.4) 
Then the initial wave functions are 

|V'(o)) = {o,...,|v^,),...,o} 

~ X/m ^nifmn'Pm (3-5) 

When the vortex is moving, the impurity sites within the core are modified and the 
Hamiltonian H. becomes time dependent. We investigate the time evolution of \ip{t)) which 
is expressed in the basis g of the instantaneous Hamiltonian 7i(t) : 

n{t)gr,{t) =E,,{t)grXt) (3.6) 

IV^n(t)) =T.m9mn4>m (3.7) 
fn — ^rn^mn(t)gm (3-8) 
(^mni,^) 9m ' fn (3-9) 

As the vortex moves at a constant velocity v, the time-dependent spectrum shows dif- 
fusive time-evolution. This is the characteristic behavior which suggests that consecutive 
Landau-Zener transitions occur incoherently. The energy transfer takes place at the points 
called avoided crossings, where the separation between levels is much less than the mean 
level spacing A. We assume that the quantum interference term in the transition probability 
between the consecutive levels can be ignored. The full quantum coherent time-evolution, 
however, leads to localization in energy space. The condition for this assumption is that 
many Landau-Zener transitions contribute to the diffusion, so we may take enough statistical 



average. [Tj 



Thus the occupation probability |?/'„(t)p spreads diffusively as time. Following ref. JT7 
we define the probability distribution Ili{t) as 

U,it) = J2ig,-f,)ik-ty (3.10) 

k 

Ili{t) is averaged over samples (over different configurations of initial impurity sites), and/or 
over many initial states i. The diffusive behavior of < n(t) > for small t is expected to be 
linearly dependent on time. 
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(n(t)) = Rt. 



(3.11) 



Solving the time dependent Schrodinger equation numerically, we investigate the dependence 
of the diffusion constant R on the rate of change of the external perturbation, i.e. on vortex 
velocity v. 

The diffusive behavior is also dependent on the smoothed density of states p. In this 
study we assume that p is constant as p = where A is the mean level spacing. There 
are two different regimes which are characterized by the small v or the large v respectively. 
The calculation of the diffusion constant R for the levels far from the Fermi energy reveals 
the nonlinear behavior of the energy dissipation of quasi-particles inside the vortex core in 
the low V regime. 

In fig. 3 we plot the result for the clean phase iVj = 3 and the dirty phase Ni = 200 
for various values oi v = 0.001,0.002,0.005,0.01,0.02,0.05,0.1 and 0.2. The initial state 
is Ei = £^35 and Ei = Eio, respectively near the hump of C(0) plot for = 1 case which 



is shown by the insets. The horizontal axis is the rescaled vortex velocity k = yC{0)v. 
In the calculation of the numerical solution of eq.( p.3[ ), the discrete time step is taken as 
At = 0.002 and the sample averaging has been done over 32 different initial configurations 
of the impurities. The gradient of the dotted line which is the least square fit is 0.67. We 
find that R (x k"' where a ~ 2/3 for small v from the fit. Since the vortex velocity is written 
as V = E/B, this result is consistent with the behavior of the dissipative part of the current 
density which is obtained in ref. as 

. _ qpnimp eI^^ f E 

where Oq is the distance from the "dangerous region" to the vortex center and 0o = vr/e is 
the flux quantum. This formula is valid for v <C vf{A/Ep)'^ and in the super-clean limit 
Ni = 1 case. Since C(0) is suppressed with the decrease of uqt (or increase of Ni), the 
random matrix prediction is irrelevant as well as in the large Ni case for this low v regime 
as expected. We also note that for larger values of v, we obtain larger values of a ~ 0.9 in 
all dirty cases of A^, = 50, 100 and 200 for which we have done the calculation. 

In the meanwhile for the large v regime the random matrix theory predicts the diffu- 
sion constant R for the perturbed Hamiltonian whose matrix elements are taken from the 
Gaussian unitary ensemble (this is the relevant for the current study) as 

R = 7rC{0)X'^ (3.13) 

where C(0) is given by eq.( p.l6 ) and X = v. The vortex damping coefficient t] is defined as 

dW/dt = r]v'^ (3.14) 

where l.h.s is the rate of the energy dissipation inside the core and t] is related to the 
dissipative conductivity as axx = rjec/nB. From the quasi-classical result in the ujqt <^ 1 
case, the damping coefficient t] is given by [|l| 

V = -Tm2D^^QT (3.15) 

where n2D is the density of electrons. Thus the random matrix theory predicts [ffl 
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C(0) = n2Di^^0T). 



(3.16) 



In fig. 4 we plot C(0) v.s. Uqt for various v values. In the large v case the relation C(0) oc uqt 
is explicitely observed. 

IV. DISCUSSION 

In this study we have investigated numerically for the universal behavior of the distri- 
bution of the rescaled level velocities for the excitation spectra within the moving vortex 
core in layered superconductors. In the dirty phase where Uqt <C 1 we obtain the Gaussian 
distribution consistent with the result of the theoretical study for the nonlinear a model 
in the various disordered quantum systems. Contrary to this, the non-Gaussian behavior 
which is considered as the characteristic property which manifests the localization of wave 
functions inside the core is observed within the clean phase UqT ~ 1. 

As to the energy dissipation within the core due to the Landau-Zener transition between 
levels which occurs at the points called "avoided crossings", the nonlinear dependence of 
the axx on the external perturbation (the vortex velocity) is investigated. This is explicitly 
related to the nonlinear I-V characteristics in the mixed state in low temperature region. 

The dissipation properties are classified into three different regimes characterized by two 
parameters, the number of impurities inside the core iVj (or equivalently the inverse of ujqt), 
and the vortex velocity v. Case i)Large v and large Ni regime; in this regime the prediction 
of the RMT is valid and at the same time the behavior becomes consistent with those 
obtained from the approach within the Kubo-Greenwood formula. Case ii) Small v and large 
Ni regime; in this regime the nonlinear behavior of the axx is observed. Since the variance 
of the scaled level velocity C(0) is suppressed as the increase of N^, the random matrix 
theoretical result is irreverent in this regime no matter how becomes large. However, 
we consider that a crossover behavior for the dissipation properties within the core should 
exist in this regime, which has been found both theoretically and numerically in the case 
of the density of states. Thus the third regime should be the crossover regime denoted 
as iii) Large and moderately large v regime. We consider the result of large a values 
which are obtained in fig. 4 in large v case implicitly manifest this crossover behavior. The 
crossover properties in this regime from the super-clean result of Larkin-Ovchinikov to 
the dirty case where RMT result should become valid is the issue of future study. 
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FIGURES 



FIG. 1. A snap shot of the perturbed spectra inside a vortex core where the external pertur- 
bation is due to the motion of the vortex at velocity v. The number of impurities inside the core 
is Ni = 3. The discretized time step is At = 0.06 and v = 0.1. The size of the Hamiltonian matrix 
is 100 X 100 and 10 levels above and below the Fermi level Ei = are shown. 

FIG. 2. The averaged distribution of the level velocities P(ei) for the vortex velocity v = 0.002 
with various numbers of impurities inside the core. The solid line corresponds to the Gaussian 
distribution. The inset shows the dispersion of level velocities C(0). 

FIG. 3. Plot of dissipation constant R versus the normalized vortex velocity k for Ni = 3 and 
Ni = 200 cases. The result for Ni = 200 is shifted downwards in order to see the coincidence of 
the slope of the fitting lines in lower k regime. The time step is taken as At = 0.002. The gradient 
of the dotted line is a = 0.67 (small k) and the dashed line a = 0.89 (large k). 

FIG. 4. Graph of log-log plot of C(0) vs wqt (= l/Ni) for various values of vortex velocity v 
for At = 1. 
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